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Abstract: We propose a method for combining QCD matrix elements and parton 
showers in Monte Carlo simulations of hadronic final states in e~^e~ annihilation. The 
matrix element and parton shower domains are separated at some value yi^i of the jet 
resolution, defined according to the /cT-clustering algorithm. The matrix elements 
are modified by Sudakov form factors and the parton showers are subjected to a veto 
procedure to cancel dependence on yi^i to next-to-leading logarithmic accuracy. The 
method provides a leading-order description of hard multi-jet configurations together 
with jet fragmentation, while avoiding the most serious problems of double counting. 
We present first results of an approximate implementation using the event generator 
APACIC++. 
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1. Introduction 

The Monte Carlo simulation of multi-jet hadronic final states is a challenging problem 
that has great practical importance in the search for new physics processes at present 
and future colliders. For example, the accurate simulation of 4-jet backgrounds was 
a central issue in the search for the Higgs boson at LEP2, and multi-jets will be a 
key ingredient in signatures of supersymmetry at the LHC. 

Two extreme approaches to simulating multi-jets can be formulated as follows. 
One can use the corresponding matrix elements, which are available at leading, or in 
a few cases next-to-leading, order in ctg, with bare partons representing jets. Alter- 
natively one can use the parton model to generate the simplest possible final state 
(e.g. e^e~ qq) and produce additional jets by parton showering. 

In the matrix-element approach, a full simulation of the final state is impossible 
unless one adds a model for the conversion of the produced partons into hadrons. Any 
realistic model will include parton showering, and then one has the problem of extra 
jet production during showering and potential double counting of jet configurations. 
On the other hand the pure parton shower approach gives a poor simulation of 
configurations with several widely separated jets. 
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The interfacing of matrix-element and parton-shower event generators is a topic 
of great current interest [1-4]. For earlier work on combining these approaches see 
Refs. [5-12]. Here we suggest a method in which the domains of applicability of 
matrix elements and parton showers are clearly separated at a given value yini of 
the jet resolution variable ycut, defined according to the /cT-algorithm [13, 14] for jet 
clustering (sometimes called the Durham algorithm). Recall that two objects i and 
j are resolved according to the /c^-algorithm if 

Vij = 2 mm{El E]}{1 - cos %)/g' > yeut (1-1) 

where E^j are the energies of the objects, 9ij is the angle between their momenta and 
Q is the overall energy scale (the cm. energy in e+e~ annihilation). Two objects that 
are not resolved are clustered by combining their four-momenta as — Pi + Pj. 

The method we propose has the following features: At ycnt > Vim multi-jet 
cross sections and distributions are given by matrix elements modified by Sudakov 
form factors. At ?/cut < Vim they are given by parton showers subjected to a 'veto' 
procedure, which cancels the t/ini dependence of the modified matrix elements to 
next-to-leading logarithmic (NLL) accuracy. 

Note that we do not attempt to give a complete description of any configuration 
to next-to- leading order (NLO) in ccg, which is why we refer to "combined" rather 
than "matched" matrix elements and showers. Procedures to combine parton show- 
ers with the matrix element corrections due to the first (i.e. at the first relative order 
in ttg) hard multi-jet configuration were considered in Refs. [5-7]. Such procedures 
might be improved by including first-order virtual corrections (see Refs. [9-12]). For 
the present, our main objective is to describe any hard multi-jet configuration to 
leading order, i.e. 0{d'^^^) for n jets in e^e^ annihilation, together with jet fragmen- 
tation to NLL accuracy, while avoiding major problems of double counting and/or 
missed phase-space regions. 

In the present paper we consider the case of e'^e~ annihilation only. In Sect. 2 
we recall the NLL expressions for e+e" jet rates, and show how they can be used 
to develop a systematic procedure for improving the tree-level predictions of multi- 
parton configurations above some jet resolution y^ja. Then in Sect. 3 we show how to 
combine these modified matrix-element configurations with parton showers, in such a 
way that dependence on y-^^i is cancelled to NLL precision. In Sect. 4 we show results 
of an approximate Monte Carlo implementation of the above scheme, and finally in 
Sect. 5 we present brief comments and conclusions. 

2. Modified matrix elements 

2.1 NLL jet rates and Sudakov factors 

The exclusive e^e~ n-jet fractions at cm. energy Q and /cr-resolution 

yini = Ql/Q'' (2.1) 
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are given to NLL accuracy^ for n = 2, 3, 4 by [14] 

R2{Q^,Q)^[A,{Q^,Q)f , (2.2) 
RsiQ^,Q) =2[A,iQ,,Q)f ("^ dqV,{q,Q)A,{Q,,q) , (2.3) 
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R4iQi,Q)^2[A,{Q,,Q)] 

("^ dq r,(g, Q)Ag{Q,, q) dq' r,(g', Q)Ag{Q^, q') 

JQi JQi 

+ ["^ dqT,{q,Q)A,{Q,,q) [' dq'T,{q',q)A,{Q^,q') 
JQi JQi 

+ ["^ dqT,{q,Q)Ag{Q,,q) P dq' V f{q')Af{Q,,q')\ (2.4) 
JQi JQi J 

where are q ^ qg, g gg and g ^ qq branching probabilities 

r,(5,Q) = — — (^ln--jj (2.5) 

rK,) = $^f^, (2.7) 

Ci? = (iV^ — l)/2iVc and = for N^. cofours, Nf is the number of active flavours, 
and Aq^g are the quark and gluon Sudakov form factors 

\iQi,Q) = exp j2dqTq{q,Q)^ (2.8) 

AgiQu Q) = exp ^- J^^ dq [Tg{q, Q) + Tf{q)]^ (2.9) 

with 

A/(gi, Q) = [A,(Qi, g)]' /A,(gi, q) . (2.10) 

The QCD running coupling as{q) is defined in the MS renormalization scheme. Part 
of the contributions beyond NLL order can be included in the calculation by using 
the definition of as{q) in the bremsstrahlung scheme of Ref. [15]. 

The Sudakov form factors Aj(gi, Q) for i = q,g represent the probability^ for a 
quark or gluon to evolve from scale Q to scale Qi without any branching (resolvable at 
scale Qi). Thus R2 is simply the probability that the produced quark and antiquark 

^By NLL accuracy, we mean that the leading and next-to-leading logarithmic contributions 
ag In^" Q/Qi and In'^"^^ Q/Qi are included in the expressions for Rn{Qi,Q)- 

^The NLL approximate expressions in Eqs. (2.5,2.6) can lead to Aj > 1. In that case one should 
replace A, > 1 by 1. 
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Figure 1: Branching structure of three-jet final state. 



both evolve without branching. More generally, the probability for a parton of type 
i to evolve from scale Q to q > Qi without branching (resolvable at scale Qi) is 
MQ,,Q)/A,{Q,,q). 

In the expression (2.3) for R3, a gluon jet is resolved at scale q where 

min{?/gp, Vqg} = gV<5^ ■ (2-11) 

Recall that in coherent parton branching the evolution variable is the emission angle 
[16] and the corresponding scale is the parton energy times the angle [17]. In the 
contribution depicted in Fig. 1, the energy and angular regions of the phase space 
that dominate at NLL order are Q Eq Eq > Eg and 1 ~ > 9qg. The quark 
evolves from scale EqOqg ~ Q to Qi without branching, while the antiquark evolves 
from Eq9gg ~ Q to g ~ Eq9qg and then branches. The resulting antiquark evolves 
from q to Qi, while the gluon evolves from q ~ EgOqg to Qi, both without branching. 
Thus the overall NLL probability is 

A,(gi, Q)^M^^M Q)\{Qi, q)\{Qi, q) = r,(g, Q) FqqgiQ,, Q; q) (2.12) 

^qWl, Q) 

where the 'Sudakov factor' Fgqg is 

FqqgiQi, Q', q) = [AgiQ,, Q)f Ag{Q,, q) . (2.13) 

Taken together with the contribution in which the quark branches instead of the 
antiquark, this gives Eq. (2.3) after integration over Qi < q < Q. 

For four or more jets, there are several branching configurations with different 
colour factors. The first term in the curly bracket of Eq. (2.4) comes from Abelian 
(QED-hke) contributions such as Fig. 2, with associated probability 

^^%^rg{q, Q)Aq{Q„ q)Ag{Q„ q)^f^^rq{q', Q)A,(Q„ q')Ag{Q„ q') 

= Fgiq, Q) Tg{q', Q) Fgggg{Q,, Q ?') (2.14) 

where the Sudakov factor is now 

FqgggiQi, Q] q, ?') = [Ag{Qi, Q)f A,(Qi, q)Ag{Qi, q') . (2.15) 
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The second term in the curly bracket of Eq. (2.4) comes from contributions with 
a,q ^ qg branching at scale q followed hy g ^ gg at scale q' (Fig. 3). The probability 
of this is 

Q) lf(Q;,^) ^^(g' gl fj(Q^';^.) r,(g-, q)A,{Q„ q')A,{Q„ q') 

= r,(g, Q) Vg{q\ q) F.-.M^, Q; q, q') (2.16) 

where the factor Fqqgg is the same as that given in Eq. (2.15). 

The final term in Eq. (2.4) corresponds to diagrams like Fig. 3 except that the 
branching at q' is g ^ qq instead of g ^ gg. The factor of Tg{q',q) is replaced by 
Tf{q') given by Eq. (2.7), and Ag{Qi,q') becomes Af{Qi,q') given by Eq. (2.10). 
Thus the Sudakov factor becomes 

FggM, Q; q, q') = [A,(Qi, Q)f Ag{Q,, q)Af{Q,, q') . (2.17) 

We see that in general the overall Sudakov factor depends on the nodal values of 
the /cy-scale q, q', ■ ■ ■ at which branching occurs, and on the types of partons involved. 
There is an overall factor of [Aq{Qi,Q)]^ coming from qq production at scale Q, a 
factor of Ag{Qi, q) when a gluon is emitted at scale q, and a factor Af{Qi, q) when a 
gluon branches to quark- antiquark at scale q. Although we have explicitly discussed 
only the n — 2,3,A jet rates, this structure of the Sudakov factor is valid for any n, 
as can be derived from the generating function given in Ref . [14] . 
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2.2 Matrix element improvement 

We can improve the description of the 3-jet distribution throughout the region 
Vqq > yqgjUqg > Uini usiug the fuU tree-lcvcl matrix element squared in 
place of the NLL branching probability Tg{q, Q) in Eq. (2.12). More precisely, we gen- 
erate qqg momentum configurations according to the matrix element squared, with 
resolution cutoff t/ini = Q\IQ^ ■, and then weight each configuration by the Sudakov 
factor Fqqg{Qi,Q;q) in Eq. (2.13), where q is given by Eq. (2.11). For consistency 
with Eqs. (2.5)-(2.7), we should also use q as the argument of the running couphng 
in the matrix element squared. 

Similarly in the four-jet case of Eq. (2.14) the product rq{q,Q)rq{q' ,Q) is an 
approximation to the full matrix element squared lAiqqggl"^ in the kinematic region 
where y,,,, and ygg/ are the smallest interparton separations. Thus it is legitimate in 
NLL approximation to replace it by lAiggggl"^ in that region. The remaining factor 
Fqqgg{Qi, Q; q, q') in Eq. (2.14) is the extra Sudakov weight to be applied. 

In general, we obtain an improved description of the jet rates and distributions, 
above the resolution value |/ini, by choosing the parton configurations according to 
the tree-level matrix elements squared and then weighting them by a product of 
Sudakov form factors. The arguments of the form factors and the running coupling 
are given by the nodal values of the /cr-resolution in the branching process, estimated 
by applying the /cT-clustering algorithm to the parton configuration. 



2.3 General procedure 

The proposed procedure for generating e^e~ — > n-jet configurations at cm. energy 
Q and jet resolution yini is thus as follows: 



1. Select the jet multiplicity n and parton identities i with probability 

(0) 

P"'\n,^) = -^ (2.18) 

where cr^'^] is the tree-level e^e~ — > n-jet cross section at resolution y^-as, — 
Q\/Q'^, calculated using a fixed value as{Qi) for the strong coupling. The 
label i is to distinguish different parton identities with the same multiplicity, 
e.g. i = qqgg or qqqq for n = 4. is the largest jet multiplicity for which the 
calculation can realistically be performed (A^ ~ 6 currently). Errors will then 
be of relative order a^~^. Ideally, one should check that any given result is 
insensitive to N. 

2. Distribute the jet momenta according to the corresponding n-parton matrix 
elements squared |Aln,iP, again using fixed as{Qi)- 
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3. Use the /o'x-clustcring algorithm to determine the resolution values 1/2 = 1 > 
y3 > . . . > Hn > Umi at which 2,3, ... ,n jets are resolved. These give the nodal 
values of qj = Q^Jyj for a tree diagram that specifies the /cT-clustering sequence 
for that configuration. 

4. Apply a coupling-constant weight of ol^{(Iz)ol^{(1^ ■ ■ •Qis(?n)/[cts(Qi)]"^~^ < 1- 

5. For each internal line of type i from a node at scale to the next node at 
qk < qj, apply a Sudakov weight factor Ai{Qi,qj)/Ai{Qi,qk) < 1. For an 
external line from a node at scale qj, the weight factor is Aj(Qi,gj). This 
procedure gives the overall Sudakov factors Fi{Qi, Q;qz, . . . , qn) of Sect. 2.1. 

6. Accept the configuration if the product of the coupling-constant weight and the 
Sudakov factor is greater than a random number TZ e [0, 1] times^ [Aq((5i, Q)Y- 
Otherwise, return to step 1. 

Note that the weight assignment is a fully gauge-invariant procedure relying only 
on the types (quark or gluon) and momenta of the final-state partons. The weight 
factor is actually independent of the detailed structure of the clustering tree and is 
the same as that for the Abelian (QED-like) graph with the same nodal scale values: 
see, for example, Eqs. (2.14) and (2.16). 

An advantage of the above procedure is that it adjusts the jet multiplicity distri- 
bution to include the Sudakov and coupling-constant weights, without the need for 
separate numerical integrations. To prove this, note that the probability of accepting 
an (n,i)-parton final state, once selected, is Pn,i = (Jn,i/<^n\^ where j includes the 
weight factors. The overall probability P{n,i) of selecting an (n, i)-parton state is 
the probability of rejecting any state any number of times before finally accepting 
the (n,i) state. Thus 



p{n,z) = y: 

m=0 



-k=N 1 



p{0) 



^k=N p(0) ^k=N . ' 

l^k,j ^k,jPk,j ^k,j '^k,3 

as required. 

In the clustering step 3, attempted clustering of partons will sometimes be 
'wrong': for example, a qqg final state may be clustered first as {qq)g. The nodal 
value for the {qq) clustering is irrelevant to NLL accuracy since there is no associated 
soft or coUinear enhancement. Hence the optimal procedure is to forbid such a clus- 
tering and continue until either {qg) or {qg) is clustered. In more complicated cases. 



^Multiplying by [Ag((5i, Q)]^ increases the efficiency of the procedure, since this constant factor 
is always present. 
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e.g. qqqq., the clustering {qq) is allowed but {qq) and (gg) should always be forbidden. 
This is simply achieved by moving to the pair of objects with the next-higher value 
of ijij whenever the lowest value belongs to a forbidden combination. 



3. Vetoed Part on Showers 

3.1 Angulair ordering and veto procedure 

Having generated multi-jet distributions above the resolution value yini according to 
matrix elements modified by form factors, it remains to generate distributions at 
lower values of Ucut by means of parton showers. This should be done in such a way 
that the dominant (LL and NLL) dependence on the arbitrary parameter yini cancels. 
Any residual dependence on 7/ini could be exploited for tuning less singular terms to 
obtain optimal agreement with data. 

Note that |/ini must set an upper limit on interparton separations Hij generated 
in the showers. Otherwise the exclusive jet rates at resolution yini could be changed 
by showering. At first sight, this might suggest that we should evolve the showers 
from the scale Qi — Q^JyhA instead of Q. However, this would correspond to using 
transverse momentum rather than angle as the evolution variable, and therefore it 
would not lead to cancellation of the dependence on Inyini. 

Consider, for example, the 2-jet rate at resolution = Q^/Q"^ < yim- If we start 
from R2 at scale Qi and then evolve from Qi to Qq, we obtain a 2-jet rate of 

[A,iQ,,Q)A,iQo,Qi)f (3.1) 

instead of the correct result 

R2iQo,Q)^[\iQo,Q)f . (3.2) 

This is because, although the y.^ values in the showers arc limited by ?/ini, the angular 
regions in which they evolve should still correspond to scale (energy times angle) Q 
rather than Qi. Consequently we should allow the showers to evolve from scale Q 
but veto any branching with transverse momentum q > Qi, i.e. the selected parton 
branching is forbidden but that parton has its scale reset to the current value as an 
upper limit for subsequent branching. 

The 2-jet rate at any scale Qq < Qi is now given by the sum of probabilities of 
0, 1, 2, . . . vetoed branchings (represented by crosses in Fig. 4) and no actual resolved 
branchings. The sum of these probabilities for the quark line is 

A,(gi,g)A,(go,g)(l+ f"^ dqr,{q,Q) + f"^ dqr,{q,Q) r dq'r,{q',Q) + --\ 

( JQi JQi JQi ) 

= A,{Q,,Q)A,{Qo,Q)exp(j^^dqT,{q,Q) \ . (3.3) 
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Figure 4: Vetoed showers on two-jet contribution. 




Figure 5: Vetoed showers on contribution with two jets at scale Qi and three at scale Qo- 




Figure 6: Vetoed showers on contribution with three jets at scales Qi and Qq. 

Comparing with Eq. (2.8), we see that the series sums to 1/ Aq{Qi,Q), canceUing 
the yini dependence and giving Ag{QQ,Q). Similarly for the antiquark line, so that 
the product does indeed give Eq. (3.2). 

For the 3-jet rate at scale Qo < Qi there are two possibilities: either the event 
is a 2-jet at scale Qi and then has one branching resolved at scale Qq, or it is a 
3-jet at scale Qi and remains so at scale Qq. The first case is depicted in Fig. 5. Its 
probability is 

[ 'dqr,{q,Q)A,{Qo,q) (3.4) 
JQo 

while that of the second case (Fig. 6) is 



2[A,{Q,,Q)f 



A,(Qo,Q) 
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■7;.,r,(,,g)A,(g,,)|f^. (3.5) 



Ag(go,g) 

\{Qu Q)\ Jq. ^ ''''' -^-^v-- ^'A,iQ,, q) 



The sum is indeed |/ini-independent and equal to Rz{Qo,Q) as given in Eq. (2.3). 
Similarly for higher jet multiplicities. A general proof of the cancellation of t/ini- 
dependence to NLL accuracy is given in Sect. 3.3. 

3.2 Initial conditions for showers 

Notice in Eq. (3.5) that the vetoed parton shower from a gluon created in a branching 
at scale q > Qi starts at scale q rather than Q or Qi. On the other hand, the shower 
from the quark line starts at scale Q. In general, each vetoed shower on an external 
parton line must start at the scale value of the node at which that parton was 
'created', in order to cancel the Qi dependence of the associated Sudakov factor. In 
the case of the branching g gg, the softer of the two gluons should be regarded as 
the one 'created', the harder one being traced back to a node at a higher scale. 

The correct treatment of the branching g qq is more subtle, although less 
crucial because this branching contributes only at NLL level. The associated factor 
Af{Qi,q') in Eq. (2.4) is a correction rather than a form factor, representing the 
conversion of a gluon jet into two quark jets at scale q'. Consequently the optimal 
treatment would be as follows: for a qq pair clustered at scale q', coming from an 
internal gluon line 'created' at scale q > q', one should generate a vetoed shower from 
the gluon starting from scale q and evolving the harder gluon at each branching'' down 
to scale q', then switch to separate showers from the quark and antiquark starting at 
scale q'. If this seems unnecessarily complicated for a next-to-leading contribution, 
one may instead consider treating the quark and antiquark as being 'created' at the 
higher scale q of their parent gluon. Then the colour factor which should be Ca 
between scales q and q' is approximated by 2Cf, an error of relative order l/N^ in a 
contribution that is already non- leading with respect to Inyini. 

3.3 Proof of cancellation of y-mi dependence 

Here we make use of the generating function formahsm and results of Ref. [14] to 
prove the cancellation of yini-dependence at NLL order. Recall that the NLL jet 
fractions at /cy-resolution yini = Q\/Q^ in a quark jet initiated at scale Q are given 

by 

i^i^Hyini = gW) = ^ (^^ j UQ^.Q:u,^u,^u)i^, (3.6) 

where the quark-jet generating function (pq is [14] 



,(Qi,Q;Mg,Mg) = Mgexpj^ dqTg{q,Q) [(i)g{Qi,q;Uq,Ug) - 1] 



(3.7) 



^The softer gluon, on the other hand, is allowed to evolve down to the shower cut-off Qq. 
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0g being the corresponding gluon-jet generating function. Now we wish to generate 
the jet fractions at some lower resolution value ycut = Qq/Q^ < Vim- This is to 
be done by replacing Ui everywhere in Eq. (3.7) by a modified generating function 
<l>iiQo, Qi-i Q] Uq, Ug), representing the vetoed parton shower. To have the correct jet 
fractions at scale Qq we require that 

(piiQl, Q] 4>q, <Pg) = 4>i{Qo, Q] Uq, Ug) . (3.8) 

Consequently we must have 

(l>q(Qo,Q;Uq,Ug) = 4)q(Qo,Qi,Q;Uq,Ug)ex.pi^J^ dqTq{q,Q) [(l)g{Qo,q;Uq,Ug) - 1]| • 

(3.9) 



Hence 

^q{Qo, Ql, Q; Uq, Ug) = (j)q{Qo, Q] Uq, Ug) CXp | - dq Tgiq, Q) [(j)g{Qo, q] Uq, Ug) - 1] 



= Uqexpi^J^ dqTq{q,Q)[(j)g{Qo,q;Uq,Ug) -1]^ , (3.10) 

using Eq. (3.7) with Qi replaced by Qo for (l)q{Qo,Q;Uq,Ug). Thus the modified 
generating function (f)q{Qo,Qi,Q',Uq,Ug) differs from the full generating function 
(pq{Qo,Q',Uq,Ug) only by having Qi as the upper limit on the g-integration in place 
of Q, i.e. by having a veto, q < Qi- Note that Q remains Q in the integrand Fq, so 
this is not equivalent to an unvetoed secondary shower starting at scale Qi. Note 
also that Q is the initial scale of the quark-jet generating function in Eq. (3.7): as 
pointed out in Sect. 3.2, this is the scale value of the node at which the external 
quark is 'created'. 

A similar result holds for gluon jets. The only difference between quark and 
gluon jets concerns the treatment of the branching g ^ qq, as discussed in Sect. 3.2. 

3.4 Colour structure 

The vetoed shower from each parton evolves in the phase space for angular- ordered 
branching [18]. This depends on the colour structure of the matrix element. As 
illustrated in Fig. 7, the angular region for parton i is a cone bounded by the direction 
of parton j (and vice- versa), where i and j are colour-connected. The upper limit 
on the scale in the vetoed shower for each parton is given by the energy of that 
parton times the relevant cone angle. This prescription identifies the cone angles for 
the 'intrinsic' radiation from each parton and is correct when the matrix element 
describes parton configurations at a hard scale Qi ~ Q. 

However, in our case some of the hard partons are produced at the scale Q, 
which is much larger than the resolution scale Qi, and the relevant cone angles are 
not set directly by the final state at scale Qi. To cancel the dependence on the 
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Figure 7: Parton shower cones. 



logarithms of Q/Qi to NLL precision, the vetoed shower has to include 'interparton' 
radiation [19], i.e. soft gluons emitted at angles that are larger than the cone angles 
for 'intrinsic' radiation. In the qqg case depicted in Fig. 1, for example, the nodal 
scale is q ^ EgOqg. The vetoed shower from the antiquark has to include not only 
gluons emitted at smaller angles 9g < 9qg but also those emitted at larger angles, 
6qg < 9g < 9qq, with energies less than Qi/9g. These soft gluons emitted at large 
angles are radiated coherently by the final-state gluon and antiquark. Thus the cone 
angle for the vetoed antiquark shower is Ogq and the initial scale is Eq9qq ~ Q. 

Notice that the starting conditions for the vetoed showers are deduced from the 
application of the /c^-clustering algorithm to the parton configurations generated 
from the modified matrix elements. It is not necessary to assign a colour structure 
explicitly to the final state at scale Qi for this purpose. The relevant colour structures 
are sampled with the correct probabilities to cancel |/ini- dependence to NLL order. 
On the other hand, if a hadronization model (cluster or string) is to be applied after 
the showers, a specific colour connection structure must be provided to the model. 

If the colour structure is not unique, colour connections can be selected according 
to their relative contributions to the matrix element squared, which are well-defined 
in the limit that the number of colours A^^ is large. Corrections to the large- A^c limit 
are normally of relative order l/N^- For high parton multiplicity, when the colour 
structure is not easily computable even at large Nc, one may use the clustering scheme 
as a first approximation in assigning colour connections. This is the procedure we 
shall adopt in Sect. 4. 

4. Results 

An approximate version of the procedure described above has been implemented in 
version 1.1 of the event generator APACIC++ [20] as follows: 
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1. Cross sections for the production of 2, 3, 4, and 5 jets according to some 
l/ini are calculated at the tree-level. The tree-level cross sections are translated 
into rates via 



For each number 3, 4, and 5 of jets, the argument of is chosen to be /tjQ^, 
where the factors ^3,4,5 are adjustable parameters chosen to reproduce the mea- 
sured jet rates. Note that this determination of the jet rates is slightly different 
from the one outlined in Sect. 2.3, for simplicity and to allow extra freedom in 
fitting the measured rates. 

2. The number of partons and their flavours are now chosen according to the 
corresponding rates in Eq. (4.1). 

3. The four-momenta of the jets are generated according to the appropriate tree- 
level matrix element. 

4. The /cT-clustering algorithm is applied sequentially until only two jets remain. 
The event is accepted with probabihty equal to the weight assigned to the 
sequence of clustering, computed as described in points 4 and 5 of Sect. 2.3. 

As recommended there, the remaining two jets arc 'forced' to be a quark- 
antiquark pair. When an event is rejected, a new configuration of momenta is 
chosen, i.e. the program returns to step 3. 

5. Next the colour configuration is chosen to be identical to the topology obtained 
in the clustering step above. 

6. Finally, parton showers are generated on external hnes according to the APACIC++ 
algorithm described in Ref . [20] , except that a veto on emission with transverse 
momentum greater than Qi is applied. In APACIC++, the evolution variable is 
virtuality and angular ordering is imposed. The initial conditions on the show- 
ers appear somewhat more restrictive than those proposed in Sect. 3.2, and so 
a slight reduction in QCD radiation is expected in this approximate implemen- 
tation of the veto procedure. 

Note that within APACIC++, more options for the steps outlined above exist, 
which are described in some detail in the manual [20]. For instance, jet rates can 
be chosen according to the NLL-rates of Eqs. (2.2-2.4), in clustering to two jets the 
configuration can be rejected if the two remaining jet flavours do not correspond to 
an quark-antiquark pair, and the colour configuration of the jets can be chosen in a 
probabilistic fashion following the prescription of Ref. [6] . 

However, we find at present that the procedure above yields the best agreement 
with experimental data. It leaves a number of parameters to be tuned, namely 



"^3,4,5 (Z/ini) 




j=3 



5 



(4.1) 
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1. The value of ctg at some reference scale. We have chosen the scale of LEP 1, 
the Z-pole. For the results displayed in the Figures, as = 0.1127 was found in 
the tune of Ref. [21]. 

2. The value of the jet resolution parameter yini at which one divides the phase 
space into a region populated by the matrix elements and the region populated 
by the parton showers. The weak (beyond NLL) dependence on this parameter 
has been employed for optimizing agreement with data. In the tune, the value 
of i/ini was fixed to 7/ini — 10~^-^. 

3. The values of the three scale factors ^3 4 5. These are supposed to compensate 
to some extent for the absence of subleading corrections to jet rates at the 
parton level. The tune gave ^3,4,5 = lO"^-^^'"^ *^®'"^-'^^. 

The parameters above together with the infrared cut-off of the parton shower and 
some fragmentation parameters have been tuned recently; for more details we refer to 
[21]. In the following we display some illustrative results, comparing the performance 
of APACIC++ with the standard event generators HERWIG [22], PYTHIA [23], ARIADNE 
[24] and with data taken by the DELPHI collaboration. The parameters of HERWIG, 
PYTHIA and ARIADNE were tuned in Refs. [25], [26] and [27], respectively. 

In Fig. 8 we depict the differential jet rates at the Z-pole as functions of the 
variable y„, which is the value of ycut a-t which an n-jet event becomes an {n — 1)- 
jet event. Clearly, all three event generators depicted here reproduce the shape of 
the distributions: deviations are on the level of at most 20% in the statistically 
significant bins. In general, APACIC++ tends to underestimate the first bins of the 
3-^2 and 4^3 distributions with an overshoot in the higher bins. This behaviour 
is somewhat reversed for the 5 — > 4 distribution. 

Integrated jet rates taken at a cm. energy of 189 GeV, defined here by the 
Cambridge algorithm [28], are displayed in Fig. 9. They demonstrate that APACIC++ 
extrapolates correctly to higher energies with all parameters fixed at the Z-pole. 

To show that the approach outlined above does indeed reproduce not only the 
correct number of jets but also the overall shape of the events, we display some event 
shapes taken at the Z-pole (Fig. 10). 

In Fig. 11 we depict some momentum spectra. Here, all the event generators tend 
to underestimate the high-momentum regions. Given the fact that the overall shapes 
of the events tend to be reproduced fairly well by the generators, one is tempted 
to conclude that this reflects a lack of particle multiplicity in the high-momentum 
regions. 

However, it should be stressed that the error bands in the right-hand plots con- 
sists of experimental errors - statistical and systematical - only. Monte-Carlo errors 
of the event generators are not included. To give some idea of the relative size of 
these errors, the numbers of events for the plots at 91 GeV are listed in Table 1. 
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DELPHI 


Ariadne 


Hcrwig 


Pythia 


APACIC++ 


350000 


2000000 


250000 


2000000 


500000 



Table 1: Number of events used to generate the plots. 



5 . Comments / Conclusions 

• Modified matrix elements plus vetoed parton showers, interfaced at some value 
t/ini of the /cT^-resolution parameter, provide a convenient way to describe simul- 
taneously the hard multi-jet and jet fragmentation regions. 

• The matrix element modifications are coupling-constant and Sudakov weights 
computed directly from the /cr-clustering sequence, which also serves to define 
the initial conditions for the parton showers. 

• Dependence on yini is cancelled to NLL accuracy by vetoing yij > yini in the 
parton showers. 

• This prescription avoids double-counting problems and missed phase-space re- 
gions. 

• In principle one needs the tree-level matrix elements |A1n,iP for |/cut > Vim at 
all values of the parton multiplicity n. In practice, if we have n < N, then yini 
must be chosen large enough for Rn>N{yim) to be negligible. 

• An approximate version of this approach (with N = 5) has been implemented 
in the event generator APACIC++ [20]. The results look promising: a rather good 
description of multi-jet observables can be achieved, and residual dependence 
on 7/ini is weak. 

• It should be possible to extend this approach to lepton-hadron and hadron- 

hadron collisions. In particular, the procedure discussed in Sects. 2 and 3 can be 
extended to deep-inelastic lepton-hadron scattering by using the corresponding 
calculation of multi-jet rates performed in Ref. [29]. 

• Extension to NLO along the lines of Refs. [9-12] may also be possible. 

Taken together, the results show sufficient agreement with data to conclude that 
this approach to combining matrix elements and parton showers is successful and 
merges the benefits of both in a rather simple way. This approach can also be used to 
introduce corrections due to the finite mass of light (with respect to the cm. energy) 
quarks, by combining the massive-quark matrix elements with the corresponding 
angular-ordered parton shower [30]. It has to be mentioned, however, that some 
significant deviations from the data remain. Therefore, additional improvements 
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such as the inclusion of NLO matrix elements seem to be necessary to achieve better 
agreement. 
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Figure 8: Differential 3 — > 2, 4 — > 3, and 5 — > 4 jet rates in the Durham algorithm at the 
Z-pole. DELPHI data (points) are compared to results (curves) of parton shower Monte 
Carlo generators. The shaded regions denote the size of the experimental errors. 
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Figure 9: Integrated 2- and 3-jet rates defined by the Cambridge algorithm at ^/s = 
189 GeV. Note that the jet rates predicted by APACIC++ are in good agreement with the 
experimental ones in the regime of the matrix elements, i.e. to the right of log^g Vini — —2.4 
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